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ABSTRACT 



We explore the consequences of an early population of intermediate mass stars 
in the 2-8 M range on cosmic chemical evolution. We discuss the implications 
of this population as it pertains to several cosmological and astrophysical observ- 
ables. For example, some very metal-poor galactic stars show large enhancements 
of carbon, typical of the C-rich ejecta of low-mass stars but not of supernovae; 
moreover, halo star carbon and oxygen abundances show wide scatter, which 
imply a wide range of star-formation and nucleosynthetic histories contributed 
to the first generations of stars. Also, recent analyses of the 4 He abundance in 
metal-poor extragalactic H II regions suggest an elevated abundance Y p ~ 0.256 
by mass, higher than the predicted result from big bang nucleosynthesis assuming 
the baryon density determined by WMAP, Y p = 0.249. Although there are large 
uncertainties in the observational determination of 4 He, this offset may suggest a 
prompt initial enrichment of 4 He in early metal-poor structures. We also discuss 
the effect of intermediate mass stars on global cosmic evolution, the reionization 
of the Universe, the density of white dwarfs, as well as SNII and SNIa rates at 
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high redshift. We also comment on the early astration of D and 7 Li. We con- 
clude that if intermediate mass stars are to be associated with Population III 
stars, their relevance is limited (primarily from observed abundance patterns) to 
low mass structures involving a limited fraction of the total baryon content of 
the Universe. 
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Introduction 



One of the outstanding questions concerning the first epoch of star formation is the 
stellar mass distribution of Population III stars. Most often, Pop III stars are assumed to 
be very massive (40 - 100 M Q ) or supermassive (> 100 M ). Tracing chemical abundance 
patterns through cosmic chemical evolution allows one to test theories employing massive 
Pop III stars. For example, the abundance patterns seen in extremely iron-poor stars do 
not support the hypothesis that the first stars had masses between 140 — 2 50 M and 
end as pair-instability supernovae (jVenkatesan et al. 1120041 ; iDaigne et al. 1 120041 ) . Yet other 
observations, such as those relating to the reionization history of the Universe, do point to 
an early population that is more massive than a standard Pop I or II distribution. Here, we 
will explore the consequences of an early population of intermediate mass (IM) (2-8 M Q ) 
stars, which dominate at a redshift z ~ 10. 

Whi le there are many theor e tical arguments poin ting; to very massive stars dominating 
Pop III (iBromm fc Larson! 12004 ; iBromm et al. 1120091 ). there exist arguments pointing to a 
first generation of stars with an initial mass function (IMF) peaked around the intermediate 
mass range. For example, the metal-free IMF predicted from opacity-limited fragmentation 
theory would peak a round 4-10 M Q with steep declines at both larger and smaller masses 
( lYoshii fc Saiolll986l ). Primordial CMB regulated-star formation also leads to the production 
of a population of early intermediate mass stars at low metallicity (Smith et al. 2009). 

On the observational side, there is abundant evidence for an early contribution from 
IM stars. Intermediate mass stars have particular effects on chemical abundance patterns. 
Unlike their massive counterparts that end as core-collapse supernovae, these stars produce 
very little in the way of heavy elements (oxygen and above), but produce significant amounts 
of carbon and/or nitrogen and above all helium. Indeed, t here is evidence tha t the number 
of carbon enhanced stars increases at low iron abundances (IRossi et al. Ill999[ ) necessitating 
a Pop HI source of carbon, possibly in the a s ymptotic giant branch (AGB) phase of IM 
stars (IFujimoto et al. Il2000t lAoki et al. 1120021 ; iLucatello et al. 1 120051 ) . To account for the 
large [C,N/Fe] ratios found in a large number of extremel y iron-poor stars in our Galaxy, it 
was argued that an IMF peaked at 4 - 10 M is needed (lAbia et al. Il200ll ). Furthermore, 
the presence of s-process elemen ts, particularly Pb at very low metalli city also points to 
AGB enrichment very early on ( lAoki et al. I l2001t ISivarani et al.l 120041 ) . As a result, an 
early population of intermediate mass stars may lead to a plausible explanation of recently 
observed carbon enriched iron-poor stars. 

Isotopic abundances of C and Mg also lend clues to the early nucleosynthetic history of 
the Galaxy. Isotopic studies of the latter show a considerable amount of dispersion at low 
metallicity ranging from low values of 25,26 Mg/ 24 Mg, consistent with pollution from core col- 
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lapse supernovae, to high values of the Mg ratios indicating the presence of AGB production 
of the heavy Mg iso topes. Core c ollapse supernovae produce almost exclusively 24 Mg, while 
the observations of lYong et al. I ( I2003af ) show enhancements (relative to predictions based 
on standard chemical evolution models) in both 25 ' 26 Mg. AGB stars were also concluded to 

Mg/ 24 Mg ratios seen in the globular cl uster NGC 67 5 2 wit h 



be the source o 
[Fe/Hl = -1.62 


f the high 25 ' 26 Mg/ : 


Yone et al. 


2003b) 


Fenner et al. 


( 


2003 


)• 



Determining the identity and nature of the first stars is a primary goal in any attempt to 
understand the physics of the first billion years of the universe. Cos mic chemical evolutionar y 
models enable one to test many of the hypotheses for Pop III stars (jCiardi fc Ferrarall2005l ). 
The bulk of the existing work in this area considers either massive stars (M > 40 M ) or 
very massive stars (M > 100 M Q ). However, given the indications discussed above for an 
early population of intermediate mass stars, it seems crucial to examine the consequences 
for cosmic chemical evolution with IM Pop III stars. 

The paper is organized as follows: In section [21 we discuss some of the expected conse- 
quences for a population of intermediate mass stars. These include a) the effective prompt 
initia l enrichment of helium, which we comp are to recent determinations of primordial he- 
lium ( llzotov fc Thuanll2010l ; lAver et al. Il2010l ) that show an excess relative to the predictions 
of big bang nucleosynthesi s (BBN); b) the reionization of the universe, usually ascribed to 
very massive Pop III stars (lCenll2003l ; lHaiman fc Holderll2003r IWyithe fc Loebll2003l ; iBromm 
20041 ); c) the evolution of the IMF; d) carbon-enhanced metal-poor stars; e) the density of 
white dwarfs; f) the rate of Type la supernovae at high redshift; and g) the evolution of N 
and Mg. 

In section [21 we introduce a relatively simple model of galactic chemical evolution to 
highlight some of the (non-cosmological) effects of the IM stars, pa rticularly on the evolu- 
tion of element abundances based on the work of iFields et al. I ( 1200 ll ). This model employs a 
bimodal star formation rate (SFR) which supplements a normal IMF and SFR, with one en- 
hanced with IM stars. In section HI we introduce a cosmic chemic a l evol ution ary model based 



heavil y on the hierarchical models developed in iDaigne et al. I (120061 ) and iRollinde et al 



( 120091 ). We present our results in section [SJ Discussion and concluding remarks are given in 
section [6j 
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2. Consequences of an early population of intermediate mass stars 



2.1. The primordial helium abundance 

Among the successes of the WMAP determination of cosmological parameters is the 
accurate determination of t he baryon density, VLnh 2 or equivalently the baryon-to-photon 
ratio 7] = n b /n 7 = lO" 10 /^; iKomatsu et~al\ hood , boioh find 



Vw 



6.19 ±0.15 



(1) 



As a consequence, it has be come possible to treat stan dard big bang nucleosynthesis (BBN) 
as a zero-parameter theory ( ICyburt et al. 1120021 120031 ) . resu lting in relatively precise predic- 
tions of the light e 
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ement abundances of D, 3 He, 4 He, 7 Li 



2004; 



Descouvemont et al. 2004; Cuoco et al. 



Coc and Vangioni 201ol ). 



t et al. 



2001 



Coc et al. 



2004! ; ISerpico et al. 112004 : 



When compared with current observational determinations of the abundances of D, 4 He, 
and 7 Li, we find varying degrees of concordance, none of which are perfect. Deuterium is 
typically credited as being responsible for the concordance between BBN and the cosmic 
microwave background determination of rj. At the current value of rj (eq. [1]), the deuterium 
abundance is predicted to be D/H = 2.52 ± 0.17 x 10~ 5 . This can be compared with the 
abundance of deuterium measured in quasar absorption systems. The weighted mean value 
of seven systems wi th reliable abundance determinations is D/H = (2.82 ± 0.21) x 10~ 5 
( iPettini et al. 1120031 . and references therein). However, the individual measurements of D/H 
shows considerable scatter (a variance of 0.53) and it is likely that systematic errors dominate 
the uncertainties. While the agreement is certainly reasonable, we do draw attention to the 
fact that the predicted abundance is somewhat lower than the observed mean and this 
(slight) discrepancy will serve as a constraint below. 

The discrepancy between t he predicted abundan ce of 7 Li and the value found in halo 
dwarf stars is well documented (ICyburt et al. 1120081 ). At the WMAP value of r] W = 6.19, 



Li/H = ( 5.12tn' 7 9) x 10 10 which is considerably higher than most observational deter- 



minations (Spite fc Spitel 19821; Spite et al. 



Bonifacio et al. 



Sbordone et al. 



& so 
2007t 



Hosford et al. 



2009 



9961: iRvan et al. 1 12000 : lAsplund et al. 1 12006 



20101). 



fl2010h finds 7 Li/H 



For example, the recent analysis of 
1.58 ± 0.31) x 10~ 10 . We comment below on the 
possible effects of astration in chemical evolution models on the lithium abundance. 

Among the light elements, helium has the most accurately predicted primordial abun- 
dance. In addition, it is relatively insensitive to the baryon density and at r] — 6.19, the 
helium mass fraction is 



Y„ = 0.2487 ±0.0002. 



(2) 
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The standard BBN results obtained in ICoc and Vangioni I (120101 ) are similar to those given 
above (Y = 0.2475 ± 0.0004, D/H = (2.68 ± 0.15) x 10" 5 , 7 Li/H = (5.14 ± 0.50) x 10~ 10 ). 

On the observational side, the det ermination of the helium abundance in extragalactic 

HII regions is fraught with difficulties ( lOlive fc Skillmanll200ll ). Over the l ast 15 years, im- 

provem ents in the analysis has led to a sy s temat ically higher 4 He abundance. Ilzotov. Thuan. fc Lipovetsky 
(119941 ) and lPeimbert. Peimbert. fc Ruizl (120001 ) introduced a 'self-consistent' method for de- 
termining the abundance of 4 He simultaneously with the determination of the physical pa- 
rameters associated with the HII region. They found Y„ = 0.240 ± 0.005 and 0.235 ± 0.002 



respectively. Using the Mont e Carlo methods described in lOlive &: Skillmanl (l200lf) and se- 
lected high quality data from Ilzotov. Thuan. fc Lipovetskyl (11994 ). lOlive fc Skillmanl (120041 ) 



found a higher value with considerably larger uncertainties, Y p = 0.249 ± 0.009. 



Further i mprovements to the analysis were aided b y ne w atomic data 



Porter. Ferland. fc MacAdam 



20071 ) used in lPeimbert. Luridiana. fc Peimbertl (120071 ) and Ilzotov et al. I (120071 ) where higher 
helium mass fractions were obtained 0.2477± 0.0029 and 0.2516 ±0.0011, respectively, based 
on the new He I emissivities and the collisional excitation of H I. The more recent work of 



Izotov & Thuanl (J2010|) found a still higher value Y p = 0.2565 ±0.0010(stat.) ±0.0050(syst.) 
Using t he same relatively small set of high quality data from Ilzotov. Thuan. &: Lipovetsky 
Jl994T ). lAver et al. I (hold ) found a similarly high value of Y p = 0.2561 ± 0.0108. The latter 
was a combined and uniform treatment of both hydrogen and helium emission lines, the 
inclusion of the new He emissivities, wavelength-dependent effects of underlying absorption, 
and neutral hydrogen collisional corrections. A combination of strong degeneracies leading 
to large uncertainties in the individual helium abundance determinations found in an unre- 
stricted Monte Carlo analysis as well as the relatively short baseline used in the regression 
with respect to O/H, the large uncertainty of 0.0108 makes certain that this determination 
is consistent with the BBN prediction at the WMAP value of r). Given, the short baseline 
of this data set, it may be argued that a weighted mean of the helium results is warranted 
rather than the commonly used reg ression of Y vs O/H as in the case of all results quoted 
above. The weighted mean found in lAver et al. I (120101 ) is 



Y p = 0.2566 ± 0.0028, 



(3) 



and one could be tempted to claim a slight discrepancy between the predicted and derived 
4 He abundance. 

If the central value of Y p is indeed 0.256, rather than invoking non-standard BBN or a 
non-standard cosmological history, it may be possible to argue for a prompt initial enrichment 
of He by an early population of stars. While the helium data responsible for Y p is at relatively 
low metallicity, the most metal-poor HII region used in the analysis leading to Eq. ([3]) is 
about 1/20 of solar metallicity. There is, therefore, a reasonable possibility that some helium 
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was produced early in the star formation history of pre-galactic structures. Here, we consider 
specific models of chemical evolution which may be responsible for such a prompt initial 
enrichment of 4 He remaining at the same time consistent with a multitude of evolutionary 
constraints. 



2.2. Reionization 



When WMAP first reported their result for a large optical depth due to electron scatter- 
ing ( IKogut et al. II2003I ). which suggested a very early epoch of reionization of the intergalac- 
tic medium (IGM), it had an immediate impact on theories of cosmic chemical evolution. 
The reported optical depth of r e = 0.17 ± 0.04 implied a redshift of reionization, z r ~ 20. 
This put an immense amount of pressure on the models of the first s t ars, and certainly 
pointed towards very mass i ve stars as a potential solut ion 



20031 : IWvithe & Loeb 



2003 



WMAP flKomatsu et al. 
to be r e 



Bromm 2004 



Cen 2003; Haiman k Holder 



Daigne et al. Il2004j ). In subsequent analyses by 



20101 ) . the optical depth was lowered and currently it is reported 
0.088 ± 0.015, corresponding to z r = 10.5 ± 1.2. This change is significant as it 



greatly relaxes the constraints imposed on the models. First, reionization at the relatively 
late redshift of 11 allows more time for structures and stars to form, and more importantly 
at the later time, more baryons are incorporated in structures (in the hierarchical structure 
formation picture), and this decreases the ionizing flux per star needed for reionization. 

The early epoch of reionization is one of the driving forces behind the identification 
of Pop III stars with supermassive stars. However, it has been argued that the number 
of very massive Pop III st ars needed for reionization would produce metal enrichments in 
conflict with observations (jVenkatesan et al. 1120041 ; iDaigne et al. 1120041 ) . On the other hand, 
it has been shown that a 'normal' population of stars (i.e., a model with a standard SFR 
and Salpeter IMF) does n ot produce enough ionization photons so as to reionize the IGM 
(IDaigne et al. 1120041 . 120061 ). Because of the short timescales involved (to achieve reionization 
by a redshift of 11), only stars with masses M > 3M can contribute to reionizing the IGM 
by a redshift of 11. While it has been shown that massive stars (in the ran ge 40-100 M^) 



are fully ca pable of reionizing the IGM in a bimodal model of star formation (IDaigne et al. 



20041 . 120061 ). to the best of our knowledge the capabilities of IM stars have yet to explored 
in this context. We will show below that indeed IM stars are capable of reionizing the IGM 
sufficiently rapidly. 
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2.3. Evolution of the IMF 



One can make a phenomenological case from galaxy evolution for an IM star-dominated 
IMF at z ^ 2. The arguments are individually refutable but the cumulative effect is certainly 
suggestive. 



Comparison of the cosmic star formation rate with the rate of stellar mass assembly, 
inferred from integrating the luminosity function over various bands, reveals a dis- 
crepancy at z ^ 2 that increases towards higher z. The discrepancy is such that one 
apparently requires a reduced ratio of stellar mass to luminosity in star-forming galax- 
ies. This is achievable if the IMF has an excess of intermediate mass stars at these 
epochs . No universal power law fit is acceptable at z £ 0.5 according to lWilkins et al. 
(j2008al ). A similar discussion also leads to a preferr ed SFR at z ^ 3 that is ~ 3 times 



higher than the best fit to the stellar mass density (IWilkins et al. Il2008bl ). There are 
other possible explanations, including uncertainties in stellar mass measurements and 
dust extinction, but no one of these is preferred. 

• Comparison of the rate of luminosity evolution of massive early-type galaxies in clusters 
at 0.02 < z < 0.83 to the rate of their color evolution requires that the IMF is 
weighted towards more massive stars at earlier epochs. A specific interpretation is 
that the characteristic stellar mass (~ 0.3M Q today) had a value of ~ 2M at z ~ 4, 



(Ivan Dokkuml 120081 ) . 



• The evolution of the galaxy stellar mass-star formation rate relationship for star- 
forming galaxies constrains the stellar mass assembly histories of galaxies. Simulations 
including gas accretion reproduce a tight correlation but the evolution with redshift 
disagrees with the data. The simulations generally prefer a constant specific star 
formation rate M*/M*. Th e specific star formatio n rate is high at z ^ 2, remains high 
to z ~ 7, bu t is low toda y ( [Gonzalez et al. Il2010al ). To reconcile these results with the 
simulations, iDavd ( 120081 ) proposed that the stellar initial mass function (IMF) evolves 
towards more high-mass star formation at earlier epochs such that the critical stellar 
mass in the IMF varies approximately from the present as (1 + z) 2 to z ~ 2. 



The arguments listed above apply to the high redshift universe where the dominant 
contributors to both light and stellar mass are massive star-bursting galaxies. There is also 
some intriguing evidence from the local universe that differentiates the IMF in current epoch 
luminous star-forming galaxies from their lower luminosity counterparts. 

Nearby massive star-forming galaxies seem to have an IMF that is more massive star- 
dominated than lower mass galaxies. There is evidence for a non-universal stellar initial mass 
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functi on from the integrated properties of SDSS galaxies, according to lHoversten fc Glazebrook 



( 120081 ). These authors inferred the stellar initial mass function from the integrated light prop- 
erties of galaxies by comparing galaxy H-alpha equivalent widths (EW) with a broadband 
color index to constrain the IMF. They find that low luminosity galaxies have much lower 
EW(H) than expected for their colors and argue that this is due to the IMF in low luminosity 
galaxies having fewer massive stars, either by steeper slope or lower upper mass cutoff, than 
that for luminous galaxies. 

A similar result for Hl-selected gala xies emerges vi a the H -alpha and the far-ultraviolet 



continuum tracers of star formation. iMeurer et al. I (120091 ) find that the flux ratio H- 



alpha/FUV shows strong correlations with the surface-brightness in H-alpha and the R 
band: Low Surface Brightness (LSB) galaxies have lower ratios compared to High Surface 
Brightness galaxies. The most plausible explanation for the correlations are systematic vari- 
ations of the upper mass limit and/or slope of the IMF at the upper end. Massive galaxies 
in the local universe have an IMF that is preferentially enhanced in massive stars relative to 
low mass galaxies. 

It seems that a case can be made that galaxies which form stars rapidly, specifically 
massive galaxies at high redshift and their nearby counterparts, possess an IMF that is 
relatively top-heavy (or bottom-light). We can speculate that essentially all vigorously star- 
forming systems at high redshift may similarly possess an IMF enriched in intermediate mass 
stars. These at least are the building blocks that in our model provide the sources of the 
bulk of ionizing photons at z as high as ~ 10 and play a role in the recycling of processed 
gas. 



2.4. CEMPs and CRUMPs and the transition discriminant Arans 

It is well known that massive Pop III stars have a very specific impact on the nucle- 
osynthesis budget of the early Universe. In this context, it is useful to incorporate stellar 
halo observations, specifically from extremely metal-poor stars (EMPS, [Fe/H] < —3), ultra- 
metal-poor stars (UMPS, [Fe/H] < —4) and hyper-metal poor stars (HMPS, [Fe/H] < —5) 
to derive constraints on the mass range of massive Pop III stars. From the nucleosynthetic 
point of view, halo stars have long been used to constrain galactic chemical evolution mod- 
els, but were, until recently, quite disconnected from cosmological models. Indeed, these old 
low mass stars preserve their pristine composition on their surface. The abundance patterns 
seen in these metal-poor stars, originating presumably from Pop III stars, can yield very 
precious information about the first nucleosynthesis processes in the Universe (subsequent 
to the Big Bang), and thus gives one the possibility of exploring the Universe during the 
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period of reionization ( IFrebel et al. 1 120071 ; lFrebelll2010l ). 



On the other hand, it has been shown (jBromm fc Loebll2003l ; IFrebel et al. 1 120071 ) that 
the abundances of ionized carbon and neutral atomic oxygen a r e imp ortant for the transi- 
tion from Pop III to Pop II /I star formation. iBromm fc Loebl (120031 ) defined a transition 
discriminant, 

Arans=log(10^+0.3x lO^) (4) 

which clearly reveals the nucleosynthetic imprint of Pop III stars. This function has the 
great advantage of being able to connect aspects of cosmo logical evolution with the C/O 
abundances in low mass stars. In this context, carbon-enhanced metal-poor stars or carbon- 
rich ultra-metal-poor stars (CEMPS and CRUMPS) seem to indicate an overabundance of 
carbon (and oxygen) at early stages of evolution. CEMPS represent about 20 percent of 
the metal poor stars. Recent studies considering multi-zone chemical evolution models show 
that only models with increased C yields or an initial to p-heavy IMF in the IM star mass 
range provide good fits of observations in the Milky Way (jMattssonll2010l ). Note that these 
models also predict a nitrogen enhancement which will be considered below. Therefore, from 
the nucleosynthesis viewpoint it is interesting to consider IM stars as Pop III candidates, 
since they are producers of both C and He. 

Finally, observations of metal absorption features in the intergalactic medium (IGM) 
give important constraints on models of the formation and ev olution of the earliest structures. 
Recent measurements of C in the IGM (IKramer et al. Il2010l ) show an increase of this element 
at about z = 6. This seems to clearly indicate an increase in the star formation rate at this 
redshift able to produce carbon. Again, IM stars are good candidates for explaining this 
trend. Note that we also predict that IGM material could have a high C/O ratio; indeed, we 
can obtain insight into this ratio from observations, though presently, it is difficult to obtain 
reliable C/O ratios from ionized abundance data (namely CIV and OVI). 



2.5. Element Abundances: N and Mg 



In the course of standard stellar evolution, nitrogen is produced through He burning. 
Primary nitrogen is also synthesized in massive stars in the H burning layer via fresh carbon 
coming from the He burning core. In intermediate-mass stars, nitrogen is produced together 
with carbon. Hot bottom-burning plays a dominant role for the synthesis of N and the 
amount of nitrogen produced depends on the efficiency and the duration of nuclear burning. 
At lower metallicity, the production of nitrogen is favored over carbon. There are indeed 
systems at low metall i city, such as 47 Tuc and M71, which show large nitrogen enhancements 
(IHarbeck et al. 1120031 ; iBriley et al. 112004 ). Abundance patterns found in these systems make 
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lt difficult to argue that the nitrogen enhancements are due to int ernal contamination, an d 
lend credence to the hypothesis of an early IM stellar population (lAshenfelter et al. 112004 ). 



The evolution of the nitrogen abundance is very sensi tive to specific parameter choices 
in the chemical evolution model (lAshenfelter et al. II2004T). Inclu ded among these are the 
uncertainties stemming from the yields ( jKarakas fc Latanzzioll2007l ) particularly in the upper 
end of the IM mass range (> 3 M ). Unfortunately, this element is difficult to observe. 
Indeed, the best measurements come from the CN BX electronic transition band which is 
not generally detected. In this context, it is difficult to confront the calculated evolution of 
this element with available observations. However, it is interesting to note that N is observed 
in the few CEMPS as shown below. 

Magnesium is produced predominantly in Type II supernovae where it is formed in the 
carbon and neon burning shells in massive stars. The dominant isotope formed in SNII is 
24 Mg. The neutron-rich isotopes 25 Mg and 26 Mg are produced in the outer carbon layer 
through a capture on neon. The yield of the heavier isotopes scales with the metallicity in 
the carbon layer and as a result, very little of these isotopes are produced at low metallic- 
ity. In contrast, significant amounts of 25,26 Mg are produced during hot bottom-burning i n 
the asymptotic giant branch phase in IM stars (jBoothroyd. Sackmann. fc Wasserburglll995l ). 
These stars are hot enough for efficient proton capture processes on Mg leading to Al (which 
decays to the heavier Mg isotopes). The neutron-rich isotopes are also produced during 
thermal pulses of the helium burning shell. Here, a capture on 22 Ne (which is produced from 
a ca pture on 14 N) lead to both 2 5 Mg and 26 Mg. The lat ter process is very important for 3 
M (jKarakas fc Latanzzioll2003l ) . Several observations (jShetronelll996l ; lYong et al. H2003al Jbl) 
indicate enhancements of the neutron-rich isotopes in low metallicity stars which could ne- 



cessitate the presence of an e arly population of IM stars (jAlibes et al. 1 12001c iFenner et al. 
20031 : lAshenfelter et al. lbo04h . 



Finally, it is interestin g; to note that Mg has been observed in ultra-faint dwarfs galaxies 
( IFrebel and Bromm Il2010l ) which are representative of the nucleosynthetic processes in the 
first galaxies. Specifically, in spite of the uncertainties related to the Mg data in these 
objects, we also note that [Mg/Fe] ratio is somew hat higher (about 1) in the se structures 
compared to the local data (about 0.3) in figure 1 of IFrebel and Bromm I ( I2010I ). Once again 
this could be a signature of the presence of IM stars in these primitive galaxies. 
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2.6. White Dwarfs 



The best current limits on the halo old white dwarf population come from a combination 
of interpreting the M ACHO and EROS micro lensing events with improved white dwarf 
atmosphere modeling ([Torres et al. II2010I ). The upper limit, including white dwarfs with 



hydrogen-deficient atmospheres, allows an increased microlensing optical depth. The upper 
limit to the contribution of halo white dwarfs is comparable to that of the thick disk, and 
accounts for at most 10% of the halo dark matter between the MWG and the LMC The 
MACHO upper limit , which has n o t bee n affected by subsequent microlensing studies, is 
approximately 20% (lAlcock et al. I (120001 ) ) . For comparison, our IM Pop III component 
gives a white dwarf mass fraction very comparable to the standard one ; it differs only by 
15 % (see Figure [3] below) . 



2.7. Supernova rates 



Among the important constraints imposed on the models of cosmic chemical evolution 
are the rates of supernovae (SN). These are intimately linked to the choice of an IMF and 
SFR and represent an independent probe of the star-forming universe. While Type II rates 
directly trace the star formation history, there is a model-dependent time delay between the 
formation (and lifetime) of the progenitor star and the Type la explosion. Since existing 
data is available only for relatively low redshifts, we discuss the predicted SN rates in the 
context of three different models and give some constraints related to the SNIa rate coming 
from the IM star mode. 

Our early IM stellar population hypothesis inevitably produces many SN at high red- 
shift. The predicted SFR and SN rates are presented below. The SFR rate corrections (for 
dust, luminosity function etc) are notoriously unreliable at z £ 8 where our IM population 
(predominantly formed earlier than z ~ 10) makes a significant contribution. The SN la 
rates provide a potentially more robust discriminant. However we note here that three issues 
help reconcile our SN rate predictions with observations of SNIa at z ^ 1. Most importantly, 
the selection efficiently of SNIa, in the absence of near infrared s e lectio n, dives rapidly to 
zero at high redshift (beyond z ~ 1.5). Indeed, from lPerrett et al. I (120101 ). the extrapolation 
of the efficiency to z = 2 implies a steep drop (to only a fraction of a percent). Moreover 
the SNIa precursors are predominan tly in low mass galax ies, which contain sub- luminous 
SNIa, by as much as 0.5 magnitude (jSullivan et al. Il2010l ). Thirdly, it may be that at low 
metallicity ([Fe/H] < —1.1), the a bsence of iron lines driving optically thick winds may 
greatly inhibit the number of SNIa (INomoto et al. 1 120031 ) . Finally, note that there are large 
uncertainties related to the time delay between IM star formation and the resulting SNIa 
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explosion. These uncertainties add to the complexity of the question, since the predicted 
bulk of SNI a can be found at higher or lower high redshift, according to the value of this 



time delay (ITotani et al. 1120101 ) . 



3. Standard chemical evolution model 

There are only a few basic ingredients to a standard (non-cosmological) galactic evolu- 
tion model. These include: the specification of an IMF, usually assumed to be a power law 
of the form <p(m) ~ m~ x \ a star formation rate, which may be proportional to some power of 
the gas mass fraction or a specified function of time; a specification of stellar properties such 
as lifetimes and elemental yields. More complex models will include infall of gas from the 
IGM, and outflows from stellar induced winds. Certain observations such as the metallicity 
distribution in low mass stars may call for the prompt initial enrichment of metals which can 
be accomplished using bi-modal models of chemical evolution in which a standard mode (e.g. 
a Salpeter IMF and a SFR proportional to the gas fraction) is supplemented with a shallow 
IMF providing predominantly massive stars over a brief period of time (or metallicity). In 
this section, we will show that this approach cannot work for the prompt enrichment of 
4 He as the most massive stars would produce heavy elements very early and thus would not 
enhance the ratio of He to metals. In other words, the enrichment of helium would occur 
too late and would not explain the enhancement in low metallicity regions. 

In contrast, bi-modal models which are augmented with a population of intermediate 
mass stars (with mass between 2 and 8 M Q ), can indeed provide an enrichment of helium 
at low metallicity in simple galactic chemical evolution models. While the model presented 
in this section cannot directly be extended to a model of cosmic chemical evolution, it does 
provide some insight to the needed bimodality involving IM stars that will c arry over into 



a mor e sophisticated model based on hierarchical structure formation as in iDaigne et al 



f l2006f ) 



In iFields et al. I ( 120011 ). we constructed chemical evolution models which led to signifi- 
cant D astration at low metallicity and also had significant white dwarf production. These 
models were based on a bimodal construction where the 'massive' mode consisted of IM 
stars. Here we consider the effect of such models on 4 He production at low metallicity. 



We illustrate results with the simple closed-box model of IFields et al. I ( 120011 ). which 
neglects both infall and outflow of baryons. Baryonic cycling through stars is determined by 
the creation function C (m, t) 

dN+ = C{m, t) dm dt, (5) 



-14- 



which describes the mass distribution of new stars and the total star-formation rate *&(t) = 
J C (m, t) dm. We adopt a "bimodal" creation function 

C(m,t) = Mt)M™)+Mt)<h(™>), (6) 
each term of which represents a distinct star-formation mode. 

In the first term of eq. ([2D we encode an early burst of IM star formation; to do this we 
adopt the star formation rate 

V>i = Zb urs t— e-* /r \ (7) 
n 

with a burst timescale t\ = 100 Myr. The parameter Xburst controls the fraction of the 
baryonic mass Mb which is processed in the burst; we will show results for different values 
of this parameter, but clearly we must have a; burst < 1- Since our focus is on IM stars, the 
IMF for this mode is a power-law (Salpeter) form: <f)(m) oc m~ 2 ' 3 , for m € [2M O ,8M ]. 



The second term in eq. fl6]) describes normal star formation. Again, following iFields et al. 



( 120011 ). we adopt the classic form ip 2 = Ag(t)M gas , with A = 0.3 Gyr , and a smoothing 
factor g(t) = 1 - e -*/°- 5Gyr which ensures that star formation begins in the bursting IM 
mode, which then smoothly goes over to the standard mode. The normal-mode IMF is 
of the Salpeter form (p2{ m ) oc rn~ 2 ' 6 , but now with m G [0.1M Q , 100M Q ]. Both IMFs are 
normalized in the usual way J m <p2 (m) dm — 1- 

In Fig. [TJ we show the evolution of the 4 He mass fraction, Y as a function of the oxygen 
abundance (relative to solar oxygen). Curves are shown for different values of the burst 
mass fraction x burst i n eq. (J7|). In the case of normal star formation (i.e., when Xburst = 
and the IM mode is absent), as O increases we see He rise from its primordial value, but 
undershoots the low-metallicity data. That normal star formation fails to match the data 
traces back to nucleosynthesis patterns in high-mass stars, which dominate the production 
of O and of new He in this scenario. Specifically, for hi gh-mass stars, the eject a have 



AY/AZo ~ 1, i.e., comparable masses of new He and of O (jWoosley &;Weaverlll995l ). Thus 
at low metallicity Zq> -C Zo,© ~ 0.01, we have AY <C 0.01, whereas the data seem to require 
AFobs — 0.01. This not only explains the problem with this simple scenario, but also suggests 
the solution: invoke a stellar population which has substantially different helium and metal 
nucleosynthesis. 

Indeed, Fig. [T] shows that as we turn on the bursting IM mode via Xbm-st > 0> there is an 
abrupt rise in He at low O abundances, with a larger rise for larger Xburst- As anticipated, this 
rise now allows for the predictions to match the data, with good fits for Xbvu-st ~ 10% — 20%. 
The origin of this rise is clearly reflects the nucleosynthesis patterns of IM stars, whose 
ejecta are much more enriched in new He than in O, relative to massive stars. Specifically, 
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Fig. 1. — The 4 He mass fraction in several bimodal models of chemical evolution. Curves 
show different fractions xtmrst ( ec f E]) of baryons processed through the burst phase. The 
data point is representative of the helium data in low metallicity extragalactic HII regions. 

for IM stars at low metallicity in t he m = (4M ,8M ) we have AY/AZ Q ~ (50,300) 



(Ivan den Hoek fc Groenewegenl 119971 ). These large ratios trace back to the very small O 



production in these stars. 

The results thus far are for the simplest possible model: a single-zone closed box. Staying 
with single-zone models, even if we include infall and/or outflow the results should be similar. 
In the case of primordial infall, both the metal and new He yields will be diluted, and by 
the same factor. Thus the He-0 trends should stay the same. Similarly, in the case of 
outflows, the abundances in the gas remaining in the box are only affected if the flow ejects 
one element preferentially relative to the other. But again, given that the O and new He 
are produced in the same environments, any physically realizable winds will remove them 
together, leaving the remaining gas with He-0 trends unaffected. 

The lessons we draw from our simple model are that an increase in helium at low 
metallicity, as suggested by observations, can be accomplished at the cost of processing 



10%-20% of observable cosmic baryons through IM stars. But as shown in iFields et al. 
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( I200ll ) and discussed above, such a scenario also has substantial consequences, including 
other element abundances, white dwarfs, and Type la supernovae. To examine all of these 
issues in a more realistic manner requires a chemical evolution model in a cosmological 
context. 



4. A cosmological chemical evolution model 
4.1. Generalities 

The study of star formation in a cosmological context requires the inclusion of i) a model 
of dark matter structure formation, and ii) accretion and outflow of baryonic matter with 
respect to existing and forming structures. 



We have developed a detailed model of cosmological chemical evolution (IDaigne et al. 



20061 : iRollinde et al. Iliopgh usin g a description of nqn line ar structures, based on the standard 



Press-Schechter (PS) formalism (IPress fc Schechterlll974l ). Baryons are placed 1) within stars 



or their remnants within collapsed structures, 2) in gas within collapsed structures (the 
interstellar medium (ISM)), or 3) outside of structures (the intergalactic medium, IGM). 
The rate at which structures accrete mass is determined by a Press-Schechter distribution 
function, fp$(M, z). The model included mass (baryon) exchange between the IGM and ISM, 
and between the ISM and stellar component. In this study we assume that the minimum 
mass of dark matter halos for star- forming structures is 10 7 M . 

Once the model is specified, we can follow many astrophysical quantities such as: the 
global SFR, the optical depth, the supernova (SN) rates and the abundances of individual 
elements (Y, D, Li, Fe, C, N, O, Mg abundances) which are used to tackle specific questions 
related to early star formation and the reionization epoch. In this study, we consider two 
distinct modes of star formation: a normal mode of Pop II /I and a massive or intermediate 
mass mode of Pop III stars. 

Throughout this paper, a primordial power spectrum with a power law index n — 1 
i s assumed and we ad opt the cosmological parameters of the so-called concordance model 



flKomatsu et al. 1120101 ). i.e. fi m = 0.27, fi A = 0.73, h = 0.71 and a 8 = 0.9. 



4.2. Star formation rate scenarios 



As in the simple model described in section [3J the models considered are bimodal. Each 
model contains a normal mode with stellar masses between 0.1 M and 100 M , with a near 
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Salpeter IMF (x = 1.6) (hereafter called model 2) 



W e fit the SFR history of Pop 11/1 stars to the data compiled in iHopkins &: Beacom 
(2006) I fro m z = to 5), and to the recent measurements at high redshift by iBouwens et al. I 
( 201ol ) and Gonzalez et al. ( 2010b ). These observations place strong constraints on the 
Pop II/I SFR. The fit f or the SFR, ip(z), is based on the mathematical form proposed in 
Springel and Hernquistl (120031 ) 



V- 



aexp(b(z - z m )) 



(8) 



b + 5exp(a (z — z m )) 

The amplitude (astration rate) and the redshift of the SFR maximum are given by v and z m 
respectively, while b and b — a are related to its slope at low and high redshifts respectively. 
The normal mode is fitted using: z/ n /i = 0.3 M Q yr -1 Mpc -3 , z m ri/i = 2.6, an/i = 1-9 and 
bn/i = 1.1. The SFR of this mode peaks at z ~ 3. 

In addition to the normal mode, we add a mode for Pop III stars. Our primary interest 
here is a model for IM Pop III stars, and we choose an IM component of stars with 2-8 
M (hereafter model 3), assuming an IMF with slope 1.3. Some fraction of the IM stars 
become SNIa. The SFR parameters are the following: u ma _ = 3 M yr -1 Mpc -3 , z m m a = 10, 
onia = 1.9 and fenia = 1.1. The redshift peak is chosen so that the IM stars are born in time 
to affect low-metallicity abundances, and the normalization is chosen to allow a substantial 
effect on 4 He. 

For comparative purposes, we also consider a massive component with stars between 36 
- 100 M Q with an IMF slope of 1.6, which terminate as Type II supernovae (SNII) (hereafter 
model 1). The SFR parameters are the following: i/nib = 0.0016 M yr -1 Mpc~ 3 , z m mh 



22.8, aiTTh = 4 and bum = 3.3. For a detailed description of the model see iRollinde et al. 



( 120091 ) and lDaigne et al. I ( 120061 ). Note that both models 1 (massive) and 3 (IM) include the 
normal mode (model 2). 

In Fig. |2j we show the adopted SFR for each of the three cases considered. As one 
can see, the normal mode is chosen to fit the observations also plotted in the figure. As 
the data extend only up to z ~ 8, there is essentially no constraint on the massive mode 
which dominates at z > 20. Data at z ^ 8 is highly uncertain due to unknown systematics 
involving, among other effects, the dust corrections and adopt ed rest-frame UV lu minosity 
function, and presents essentially no conflict with our models (ILabbe et al. I (120 101 ) ). 



Fig. [2] also shows that the IM mode represents a significant processing of baryons. The 
peak of the IM star formation rate is about a factor ~ 10 higher than the peak of the normal 
star formation rate, and a factor ~ 10 3 higher than the massive Pop III star formation rate. 
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On the other hand, the IM star formation peak lasts for about a factor ~ 10 less time than 
the normal star formation rate. Thus we expect that roughly comparable amounts of baryons 
will be processed through these two modes. In contrast, a much smaller fraction of baryons 
are processes through the massive Pop III mode. One should recall, however, that in our 
model many baryons always remain the IGM and never reside in galaxies nor are processed 
though stars. 



Time (Gyr) 

13.675 1 0.5 0.1024 




Z 



Fig. 2. — The cosmic star for mation rate (SFR) a s a fu nction of redshift. The data (solid 
black points) are take n from iHopkins Beacoml (120 06). Dashed black points come from 



Bouwens et al.\ (120101 ) and from iGonzalez et al. I (j2010bl ). The blue solid line represents the 



standard SFR with a Salpeter IMF and a mass range: 0.1 < M/M <100. The dashed red 
line represents the massive Pop III stellar mode, with a Salpeter IMF and a mass range: 36 
< M/Mq < 100. The dotted black curve represents the IM SFR mode, with a Salpeter IMF 
and a mass range: 2 < M/M <8. 



Figure |3] represents the mass fraction of each component related to the total mass of 
baryons: IGM, Galaxies and stars. Note that the addition of the IM star component does 
not affect the mass fraction of stars (including remnants) at low redshift . 
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Note that the star fraction in the IM mode has an initial peak at about ~ 10% of all 
baryons. This occurs at z ~ 10. The star fraction then drops as the IM stars die off, then 
rises with normal star formation. Thus we see that in our model, about 10% of the baryons 
are processed through the IM mode. In contrast, we see that the massive Pop III mode only 
processes about ~ 10~ 6 of all baryons. 




Fig. 3. — Different mass fractions as a function of redshift. The flat solid black line (top of 
the figure) represents the total baryons in the Universe. The IGM fraction is shown by the 
dashed black curve, the galactic fraction is shown by the growing solid blue curve. The mass 
fraction of stars is plotted for the three models considered (red curves): solid line for model 
2 (normal mode), dashed line for the model 1 (with massive Pop III stars) and the dotted 
line for model 3 (with IM stars). 



4.3. Yields and lifetimes 



The lifetime s of intermediate mass stars (0.9 < M/M < 8) are taken from lMaeder and Meynet 

(119891 ) and from ISchaererl (120021 ) for more massive stars. Old halo stars with masses below 
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~ 0.9 M have a lifetime long enough to be observed today. They are assumed to inherit 
the abundances of the ISM at the time of their formation. Thus, their observed abundances 
reflect, in a complex way (due to exchanges with the IGM), the yields of all massive stars 
that have exploded earlier. 

The yields of stars depend on their mass and their metallicity, but not on their status 
(i.e. Pop II/I or Pop III). Some Pop II/I stars are massive, although in only a very small 
proportion since we use a slightly steeper than Salpeter IM F. Pop HI stars are all massive 



stars. We use the tables of yields (and remnant types) given in Ivan den Hoek fc Groenewegen 



( 119971 ) for intermediate mass stars (< 8M Q ), and the tables in IWoosley &: Weaver! (119951 ) for 
massive stars (8 < M/M < 40 ). An interpolation is made between different metallicities 
(Z=0 and Z=10~ 4 '~ 3 '~ 2 '~ 1,0 Z Q ) and we extrapolate the tabulated values beyond 40 M . 



5. Results 

Having specified the models under consideration, we can now systematically consider 
the consequences of our particular choices for Pop III stars. For most of the results which 
follow, we will compare three distinct model choices. 1) model 2 alone, ie, only Pop II/I 
stars; 2) model 3, a bimodal model including the normal mode (model 2) plus the IM mode 
of Pop III stars ; 3) model 1, a bimodal model including the normal model (model 2) plus 
the massive mode of Pop III stars. 



5.1. Reionization 

Having set the SFR in our model, we now compute the electron scattering optical depth 
for our choice of IMF. The evolution of the volume filling fraction of ionized regions is given 
by: 



dQion(z) 1 dn ion (z) 



dz 



rib dz 



- a B n h C(z)Q 2 imi (z) (1 + z)' : 



dt 



dz 



(9) 



where is the co moving density in baryons, ni on (z^ the comoving density of ionizing photons, 
an the recombination coefficient, and C(z) the clumping factor. This factor is taken from 
Greif fc Bromml ( 120061 ) and varies from a value of 2 at z < 20 to a constant value of 10 for 
z < 6. The escape fraction, / esc , is set to 0.2 for both Pop III and PopII /I . Th e number 



of ionizing photons for massive stars is calculated using the tables given in ISchaererl (120021 ). 
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Finally, the Thomson optical depth is computed as in iGreif fe Bromml ( 120061 ): 



r 



ca T n h / dz' Q ioa (z') 
Jo 



l + z 



/\3 



df 



dz' 



(10) 



where z is the redshift of emission, and <7t the Thomson scattering cross- sect ion. 



In Figure SJ we plot the integrated opti cal depth from z = to z. The red band 
represents the observed results from WMAP7 (IKomatsu et al. 1120101 ). As one might expect, 
the normal mode alone (model 2 shown as the solid blue curve) is not capable of producing 
a sufficiently large optical depth. The massive mode (model 1 shown by the red dashed 
curve) does manage to reionize the IGM, but only barely. Of course more reionization is 
possible by increasing the SFR for the massive mode, but this would lead to complications 
in elemental abundances. We see that the IM mode (model 3 - shown by the black dotted 
curve) is able to easily reionize the Universe compared to models 1 and 2. Our models allow 
an escape fraction as low as ~ 0.1 — 0.2, which provides considerable flexibility in accounting 
for reionization. While the ratio of ionizing photons from a 6M star vs that of a 2OM 
star is small, one has many additional stars in the IM scenario relative to the number in the 
massive mode. 



5.2. Nucleosynthesis evolution 



We next consider the evolution of the abundances of He, CNO, Fe and Mg as a function 
of the redshift and/or the metallicity. We also show the evolution of D and 7 Li due to early 
astration. 

We begin with the evolution of the 4 He abundance. In the left panel of Fig. [SJ we show 
the analogue of Fig. [1] for the bimodal models based on hierarchical structure formation. 
Plotted is the 4 He mass fraction, Y, vs. [O/H]. The solid blue curve corresponds to the 
model with only a Pop II/I contribution. This result is indistinguishable from that produced 
by models including the massive Pop III contribution (not shown). This can be understood 
as the bulk of the 4 He and oxygen are derived from the same stars in the two models. In 
both cases, the helium abundance begins at the BBN primordial value and remains rather 
flat until the oxygen abundance is roughy l/10th of solar. In contrast, as shown by the 
black dotted curve, in the IM Pop III model, 4 He is produced early in IM stars with little 
or no accompanied oxygen. As a result, the 4 He mass fraction begins to grow at very low 
[O/H] and for the model parameters c hosen, plateaus at value close to Y ~ 0.2 56, in good 
agreement with recent determinations (llzotov fe Thuanl |2010| ; lAver et al. 1 120101 ) . As in the 
case of the simple model discussed in section [31 this model, therefore produces an effective 
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Fig. 4. — The optical depth as a function of the r edshif t. The red range corresponds to the 
observed results from WMAP7( iKomatsu et al. I ( ]2010l )). The red dashed line corresponds 
to model 1 (with Massive Pop III stars), the blue solid line to model 2 (the normal mode) 
and the black dotted line to model 3 (with IM stars). 



prompt initial enrichment of 4 He in low metallicity galaxies. Note however, that due to the 
size of the error bar associated with the observations, we cannot exclude the possibility that 
no enrichment occurred, leaving models 1 and 2 viable. 

It is helpful to compare the helium-oxygen trend in Fig 5(a) with the simple closed box 
results in Fig. 1. The latter shows results for different fractions Xbm-st °f baryons processed 
through the IM mode. As noted above, our full hierarchical model cycles about 10% of all 
baryons through IM stars; thus Fig. 5(a) should be compared to the Xburst = 0.1 curve in 
Fig. 1. Indeed, we see that both show very similar 4 He evolution. 

In the right panel of Fig. [5j we show the corresponding evolution of the 4 He abundance 
with redshift. The dashed red curve corresponds to the massive Pop III mode. Once again, 
the standard model and the massive Pop III model show nearly identical histories. We do see 
however, a modest increase in the helium mass fraction at high redshift due to the massive 
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Fig. 5. — Left: Evolution of the helium mass fraction as a function of [O/H] (relative to the 
solar value). The solid blue line corresponds to the both models: 2 (standard model) and 1 
(with a massive Pop III mode). The do tted black curve co rresponds to model 3 ( with IM 
Pop III stars). The red point comes from lAver et al. I (j2010l ). Right: Evolution of the helium 
mass fraction as a function of redshift . The red dashed line corresponds to the Massive Pop 
III model , the blue solid line to the standard model, and the black dotted one to the IM 
Pop III model. 



Pop III mode. This is diluted by further infall as structures continue to grow. In the IM 
Pop III model, there is a significant enhancement in the helium mass fraction around z ~ 8. 
This too is slightly diluted with infall at later times. 

We also show in Figs. [6] - [7] the corresponding evolution for D and 7 Li; note that the 
vertical axes have zero offset in order to more clearly distinguish the model predictions. In 
each case, there is little difference between the standard model and the ma ssive Pop HI 
mode l. Intermediate mass stars on the other hand are known to deplete D/H ([Fields et al. 
200ll ). In this case, the degree of depletion is not severe (about 15%) but it does move the 
BBN value further away from the abundance determined in quasar absorption systems. The 
7 Li astration factor is the same as for D/H and in this case moves the BBN value closer to 
the abundance determined in halo stars. However, an abundance of ~ 4 x 10~ 10 is still very 
far from the observed plateau value between 1 and 2 xlCT 10 . 



Consider now the abundance evolution of CNOMg together with their abundance ratios. 
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Fig. 6. — Left: As in Fig. showing the evolution of the deuterium abundance as a function 
of [O/H], Right: showing the evolution of the deuterium abundance with redshift. Note that 
the vertical axis is offset from zero, in order to more clearly show the (relatively small) effect 
on D. 

Figures [8] and [9] display the abundance evolution with z. As in all of the previous figures, 
each of the following figures shows three curves corresponding to a model with no Pop III 
contribution (solid blue curve), a model with a massive Pop III component (red dashed 
curve) and a model with an IM Pop III component (black dotted curve). The evolutionary 
behaviour of these model choices is well understood. Normal mode stars produce a moderate 
abundance of heavy elements at high redshift from the IMF-suppressed massive stars in that 
mode. The same stars, nevertheless are effective at low redshift and produce the bulk of 
metals observed in the Milky Way. The Massive Pop III component produces C, O and Mg 
at high redshift, though these abundances are diluted at low redshift by infall as structures 
continue to grow. The IM stellar component produces a rather specific nucleosynthetic 
signature: essentially C and N. 

Figures fTUl and ITT1 show the evolution of the abundance ratios as a function of redshift. 
As one might expect, at lower redshifts, the IM star Pop III component produces the highest 
[O/Fe], [C/Fe] , [N/Fe]. The evolution of [Mg/Fe] corresponds to 24 Mg coming from massive 
stars and the highest value comes from the massive population III component. The ratio of 
the neutron-rich isotopes to 24 Mg remains small in the Pop III model dominated by massive 
stars, but this ratio can become large (C?(l)) in the IM Pop III model. Also shown in these 
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Fig. 7. — Left: As in Fig. showing the evolution of the lithium abundance as a function of 
[O/H] . Right: showing the evolution of the lithium abundance with redshift. 




Fig. 8. — As in Fig. |5]for Carbon (Left), and Nitrogen (Right). 

figures are the abundance ratios seen in several CEMPS. 

Surprisingly, the [C/O] ratio shown in Figure [12] is high in this model as oxygen is 
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Fig. 9. — As in Fig. |5]for Oxygen (Left) and Magnesium (Right). 

also primarily produced in more massive stars whose role is diminished in this model at 
intermediate redshifts. In contrast, the massive Pop III component reaches very high ratios 
of [C/Fe] and [O/Fe]. Only the massive Pop III model is capable of reproducing the [O/Fe] 
abundance ratios observed in these stars. Similarly, the IM model has difficulty in achieving 
[C/Fe] ratios as high as the two HE stars with [C/Fe] ~ 4 as seen in the figure, although 
the model does a reasonable job of accounting for the [C/O] ratios seen in these stars. 
The abundance patterns found in these CEMPS continues to confirm that the formation of 
massive stars at very high redshift is plausible. 

Finally we consider the evolution of the parameter D trans defined above in eq. (jl]). 
Figure [13] shows the D trans parameter as a func tion of the metal enrichment. Black circles 



come from a compilation of iFrebel et al. I (120071 ). We see that this representation allows one 



to clearly distinguish the effects of the different choices of stellar mass ranges: the normal 
mode fits the bulk of s tandard stellar data. The massive mode can explain the few CRUMPS 



( IRollinde et al. 1 120091 ) and the IM stellar mode can explain the bulk of stars very enriched 
in C between —3 < [Fe/H] < —1 range ( red box). Indeed, we can see that the black dotted 
line (IM stars) is always much higher than the blue one (normal mode). 



In particular, for [Fe/H] > —3, IM stars give D trans > — 1. In Figure 1 of IFrebel et al. 



( 120071 ) . a specific stellar population is found exactly there: C rich stars. It is very interesting 



to note that only the IM component is able to explain this part of the diagram. Indeed, 
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Fig. 10. — Left: As in Fig. |5j the evolution of the [C/Fe] ratio (relative to the solar value) as 
a function of redshift. Right: Idem for [O/Fe] ratio. Observational data (horizontal dashed 
lines ) represent measured abundances in t he following very ir on poor halo sta rs: CS 22949- 



037 



Depagne et al 1120021). He 0107-22 40 flBessel et al 1 120041 ) . He 1327-2326 (IFrebel et al 



200fih and 077-61 jpiez fc Cohenll2005l ). 



to the best of our knowledge, the yield patterns of IM stars are the only ones which can 
correctly fit the abundances seen in C stars. Thus is it tempting to interpret the metal-poor, 
C-rich stellar population as an indication of the presence of an early IM population. One 
could even push this further: of the stars in Figure [131 about ~ 10 — 20% fall in the C-rich 
red box fit by IM stars. To the extent that the data in the figure faithfully trace the ensemble 
of star-forming histories in the early universe, we would then estimate that Pop III IM stars 
form at ~ 10 — 20% of our model rate. Since our IM model has ~ 10% of baryons processed 
through this mode, the C-rich stars would in turn imply that ~ 1 — 2% of baryons participate 
in IM star formation. 

Thus we see that C-rich halo stars seem to demand some early IM star formation. 
Moreover, the large scatter in C and O abundance patterns over the full metal-poor halo 
star population demand that no single nucleosynthesis (and thus star-forming) history will 
suffice. Indeed, given that C-rich stars are outliers to the main trend, the IM population 
appears to be probably sub-dominant, acting only in low mass galactic structures at high z. 
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Fig. 11 — As in Fig. [10] for the [N/Fe] and [Mg/Fe] ratios. 
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Fig. 12.— As in Fig. Hi for the [C/O] ratio. 

5.3. Type II supernovae 



In Figure [HJ at high redshift, we show the resulting evolution of the Type II SN rate 
along with the data at relatively low redshift. The GOODS data for core collapse supernovae 
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[Fe/H] 



Fig. 13. — As in Fig. [U the evolution of the D tra ns parameter (transition discriminant D trans 



see eq. (Hj) 



Frebel et al. 



) as a function of [Fe/H] . The black circles come from the compilation of 
( 120071 ) . The red box corresponds to the region corresponding to C rich stars. 



have been placed in two bins at z = 0.4±0.2, and z = 0.7±0.2 (IDalhen et al. I (12004 120081 ) I 
Their results (which have been corrected for the effe cts of extinction) show SNII rates whic h 



are significantly higher than the local rate (at z = 0; ICappelaro et al. 111999c iLi et al. II2010I ). 



Since the timescale between star formation and the core collapse explosion is very short, 
there will be no significant contribution of Population III stars to Type II rates at any z. 
The SN Type II rate is directly related to the overall SFR and therefore to the astration rate 
v and the slope of the SFR at high redshift. 

The data at low z is well described by the models. This result is independent of our 
choice of the mass range for Pop III since massive Pop III stars will only contribute to the 
SN rate at very high redshift and our IM Pop III stars produce no SNII. However, the model 
predicts much higher SN rates at higher redshift. Indeed, the SN Type II rate peaks at 
a redshift z ~ 3 at a rate which is nearly 8 times the observed rate at low z. It will be 
interesting to see whether future data will be able to probe the SN rate at higher redshift. 
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14. — Evolution of the SNII rate as a function of redshift. The observed rates are 



taken from 


CaDDelaro et al. 


( 


1999): 


Botticella et al. 


(2008); 


Bazin et al. 



Capellaro et al. 


(2009b; 


Li et al. 



(2005): 


Dalhen et al. 


( 


2004, 


2008) 



(hoioh . The evolution is the same for 



all three models. 



These elevated rates may be detectable in spite of the expected increase in dust extinction 
due to the early production of metals. 



5.4. Type la supernovae 

Our calculation of the Type la SN rate depends on two additional assumptions beyond 
the specification of the models present in the previous section. The Type la rate will depend 
on the fraction of low and IM mass stars which end up as SNIa, as well as the time delay 
between formation and explosion. Furthermore, it is not clear that either of these quantities 
are universal constants. That is, they may vary with redshift, metallicity, or the size of the 
structure the stars are formed in. We have assumed that the SNIa rate is proportional to 
the IM star SFR (2 — 8M ). The coefficients adopted are e = 0.02 and 0.01 for the models 
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2, 1 and 3 respectively. 

In Fig. we show the SNIa rate contrasted with the observational data. As in the case 



for Type II supernovae, at high redshift, the data (also taken from GOODS; iDalhen et al. 



2004 . 120081 ) is binned into four redshift bins at z = 0.4, 0.8, 1.2, and 1.6, each with a spread 



of ±0.2. Other data references are given in the caption. We see that, in contrast to the case 
of SNII, the evolution of the SNI rate is very different for the two models with and without 
IM stars. Adding the IM stellar mode implies a high SNIa rate at high redshift. The position 
of this bump depends on the time delay of the SNIa. 





Fig. 15. — Evolution of the SN I rate a s a function of redshift. The observed rates are taken 
from ICappelaro et al. I Jl999h: Reiss 



Pain et al. 



Blanc et al I ( 120041 ): 



(2002 



Madgwick et al 



(200C 



(2003) 



Hardin et al. I (120001) : iDalhen et al 



Tonrv et al. I (12003( 1: IStroder et al 



(2004): 

(EqqJ; 



Li et al. l fcoioh : iPerrett et al. I fcoioh . The red solid curve corresponds 



to the both models 1 (normal mode) and 2 (with Massive Pop III stars). The black dotted 
curve corresponds to model 3 (with IM stars). Left: The assumed time delays for IM stars 
in model 2 and model 3 are 2.5 and 3.4 Gyr (including the lifetime of stars) respectively, and 
the fraction of the white dwarfs which become SNIa are e = 2 and 1 per cent, respectively. 
Right: The time delay for IM stars to become SNIa in model 3 is reduced to 0.5 Gyr and e 
is 2 per cent for both models. 



Due to large uncerta inties concerning the time delay, we present two cases ( IMaoz et al 



20101 : iTotani et al. Il2010l ). We hold the delay fixed at 2.5 Gyr for normal mode stars, but in 
the left panel of Fig. [151 we choose a delay of 3.4 Gyr and in the right panel we take 0.5 Gyr 
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for model 3. While the normal mode fits the data quite well, we see that the population of 
IM stars creates a significant bump in the SNIa rate at high redshift. Unless the time delay 
for these supernovae is reduced from the nominal value of 2.5 Gyr (as in the right panel of 
Fig. [15]) the SNIa rate would greatly exceed the observed rate at redshifts between 1 and 
1.5. If the delay is smaller, the bump in the SNIa rate is pushed to higher redshift. Of 
course, it may be that for Pop III stars, the fraction e of IM stars becoming SNIa is reduced 
as discussed in section 12.71 

In the upper panel Fig. [TBI the SNI/SNII ratio is plotted as a function of the redshift 
for models 1 and 3. In the lower panel of the figure both SNIa and SNII rates are plotted, 
as in [TH and [15] The IM stellar mode predicts that the bulk of SNIa occur at hight redshift, 
though this results depends on the time delay of thes e SN. However , at su ch a high redshift 



the selection efficiently drops rapidly towards zero flPerrett et al. I (120101 )) and we cannot 
exclude the existence of this SNIa component. In the future it would be interesting to have 
observational insight into low mass galaxies, which could contain sub-luminous SNIa. An 
alternative, if these sub luminous SNIa are not observed we could put strong constraints on 
the early helium abundance in the primitive structures. 

It is interesting to note the interplay between the IM constraints from halo star carbon 
abundances and from and Type la SNe. Namely, at the end of §5.21 we argued that the 
existence of C-rich metal-poor halo stars seems to require an early cosmic IM star population, 
but the relative rarity of such C-rich stars also seems to imply the IM population was sub- 
dominant, perhaps accounting for ~ 10 — 20% of early star formation. Note that if we apply 
this same ~ 10 — 20% factor to the IM predictions for Type IA SNe in Figs. [15] and [16] 
the predictions would be brought into agreement with the observations. Similarly, we could 
argue the other direction, demanding that the IM fraction be small enough to agree with 
the Type la SN observations; this would again give ~ 10 — 20% of our all-IM Pop III model, 
and thus would agree with the rough statistics suggested by the counts of C-rich metal poor 
stars. This rough concordance may be coincidental but is nonetheless both encouraging and 
intriguing. 



6. Discussion 

By definition, as the first stars formed in a metal- free environment, Pop III stars existed 
and played a role in our past history. However, the identity and mass distribution of these 
stars is largely unknown. We expect that they played an important role in re-ionizing the 
Universe at high redshift, as well as laid the chemical seeds for future generations of stars 
(population II/I). We also know that a Salpeter IMF, typical of population I with a mass 
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Fig. 16. — Evolution of the ratio of SNIa and SNII event rates as a function of redshift. (top 
panel) for model 1 and model 3. In the bottom panel, the rates plotted of SNI and SNII as 
presented in [TJ] and [151 The delay of SNIa is taken as 3.4 Gyr. 



range of 0.1 - 100 M© is not capable of producing a sufficient optical depth for reionization 
or of producing the specific abundance patterns in EMPs. Very massive stars (> 100 M ) 
also fail in producing the observed abundance patterns, though they are certainly capable 
of reionizi ng the Univers e. A t op-heavy IMF for Pop III was studied in the cosmological 



context in lDaigne et al. I ( 120061 ); iRollinde et al. I (120091 ) and while successful in many aspects, 
is not capable of producing the (albeit weak) evidence for an enhancement in 4 He at low 
metallicity over BBN predictions or the lagre carbon enhancements which are made manifest 
in the parameter -Dtrans- To the best of our knowledge, IM stars as candidates for Pop III 
have not been studied in the cosmological context. 

We have shown that an early generation of IM stars is indeed capable of providing an 
effective prompt initial enrichment of He. We have also shown, perhaps surprisingly, that 
these stars normalized in abundance to the required He enrichment, are more than capable 
of providing a sufficient number of ionizing photons to the early IGM. Also as a bonus, IM 
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stars can provide the large abundance ratios of C, N, 0, and Mg relative to Fe in EMP stars. 
However, our results indicate that one cannot simply assume a homogeneous and well-mixed 
IGM and/or ISM as smaller structures grow to large galaxies and clusters of galaxies. The 
CNO abundance that this model predicts at metallicities between — 1 < [Fe/H] < —3 far 
surpasses the abundances observed in the IGM. 

We are thus led to an interesting dilemma which is best characterized by the data 
displayed in Fig. [T3J Recall that the discriminant -D trans represents a measure of carbon 
and oxygen produced in prior generations of stars. As discussed earlier, the data show two 
distinct populations: the bulk of the data show increasing -Dtrans with [Fe/H] and is very well 
modeled by either a standard or top-heavy IMF; and a second population with very large 
-Dtrans between — 1 < [Fe/H] < —3. The latter is well explained by our Pop III IM stellar 
mode. How can one model explain both populations? 

It is perhaps too naive to expect a homogeneous model as the type we have been 
describing to explain this discrepancy. Indeed one can imagine that the impact of the 
generation of IM stars was effective only early on in the building of higher mass objects in 
the hierarchical structure formation scenario. That is, this generation was active only in the 
smallest scale structures, here taken to be typically 10 7 M Q . These structures and the low 
mass stars and remnants left behind were in some cases incorporated into the halos of larger 
objects to become galaxies. Others remain as low mass dwarfs. In this way, one can perhaps 
explain the bulk of the observations of low and increasing -Dtrans i n stars formed in larger 
scale objects involving a larger baryon fraction, while those stars with large -Dt rans would be 
explained by the impact of IM stars in the smallest structures formed. Indeed, using both 
carbon data on metal-poor stars, and Type la SN observations, we roughly estimate that a 
Pop III IM mode operates at ~ 10 — 20% of our model rate. Since our IM mode processes 
about 10% of baryons, this in turn implies that ~ 1 — 2% of baryons participate in IM star 
formation. 

In this interpretation, the prompt initial enrichment of 4 He, which is observed largely in 
dwarf irregular galaxies, was also impacted by IM stars as these objects were evidently not 
incorporated into larger scale structures. In the context of the light element abundances, it is 
interesting to note further that some of the stars which show large car bon enhancements also 



show deficiencies in 7 Li. Furthermore, there is increasing evidence (ISbordone et al. 112010 



Melendez et al. 


2010 


) for 


< -3. In 


Sbordone et al. 



Sbordone et al. I (120101 ). it is argued that the depletion (which shows considerable 
dispersion below the plateau) at low metallicity may have been due to stellar astration. The 
model discussed here (given the above interpretation) can account for most of the astration 
seen in these stars. The model cannot, however, account for the discrepancy between the 
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BB N predictio n at th e WMAP baryon density and the Spite plateau value as speculated 



in 



Piau et al. I ( 120061 ). Note that the depletion in D/H is probably not an issue, as the 
absorption systems with measured D/H are presumably larger scale structures for which our 
abundance patterns would not be expected to apply. 

Unfortunately, neither of the two Pop III models delivers entirely satisfactory results. 
While the massive mode, is capable of producing sufficiently high ratios of [C,0,Mg/Fe] to 
explain the observations of extremely iron-poor stars, the high ratios occur almost instantly 
after the first stars are born, and thus require these stars to be born at that time at very high 
redshift. This model also has difficulty in explaining the high [N/Fe] or [C/O] ratios seen. 
As a result, the massive Pop III model cannot explain the high values of -Dtrans seen in some 
stars. The model cannot account for the effective prompt initial enrichment of 4 He nor can 
it be used to account for the astration of 7 Li below the Spite plateau. On the other hand, 
the IM model for Pop III underproduces [O/Fe] and because of the large amounts of C and 
N produced cannot explain the bulk of the values of D trans observed, but the top panel data 
is explained due to the overabundances of C,N. Indeed, the overproduction of these elements 
would preclude the homogeneous and well-mixed treatment of gas in ever increasingly large 
scale structures. 

Another distinctive feature of the IM Pop III model is the prediction for the rate of 
Type la supernovae at high redshift. This model naturally predicts a large increase in the 
Type la rate if the efficiency for supernovae is constant and the time delay is not small. It is 
relevant to speculate that a small high mass tail in the IM population could lead to possible 
GRBs in the redshift range 10-20. This would fill the "GRB desert" between Pop II and the 
usual high redshift Pop III stars. 
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Thanks also to Yannick Mellier and Reynald Pain for their precious help regarding the SNIa 
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